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A quantitative description of magnons in long-range ordered quantum antiferromagnets is pre¬ 
sented which is consistent from low to high energies. It is illustrated for the generic S — 1/2 
Heisenberg model on the square lattice. The approach is based on a continuous similarity trans¬ 
formation in momentum space using the scaling dimension as truncation criterion. Evidence is 
found for significant magnon-magnon attraction inducing a Higgs resonance. The high-energy roton 
minimum in the magnon dispersion appears to be induced by strong magnon-Higgs scattering. 
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Introduction — The characterization of elementary ex¬ 
citations is of fundamental importance in condensed mat¬ 
ter physics. Especially strongly correlated systems are 
known to host a large range of fascinating collective 
behavior leading to a zoo of exotic elementary quasi¬ 
particles with integer or non-integer, i.e, fractional, quan¬ 
tum numbers. The decay of integer excitations into frac¬ 
tional ones is called fractionalization. Prominent exam¬ 
ples are the fractional quantum Hall effect with frac¬ 
tional charge excitations [ ] or the physics of magnetic 
monopoles in spin ice where spin north and south poles 
move independently [2]. Another arena of major impor¬ 
tance is high-temperature superconductivity, where the 
nature of magnetic excitations has been intensively de¬ 
bated. The key point is whether fractional spin excita¬ 
tions, spinons with S = 1/2, are present and essential [3], 
or whether a description in terms of magnons or triplons 
with integer spin S = 1 is appropriate [ ]. 

Theoretically, the physics of undoped cuprates is 
closely related to the antiferromagnetic S' = 1/2 Heisen¬ 
berg model on the two-dimensional (2D) square lattice 
[ ]. This paradigmatic unfrustrated model is known 
to display long-range Neel order with finite sublattice 
magnetization at zero temperature [6, 7]. The well- 
understood low-energy excitations are gapless magnons 
as enforced by Goldstone’s theorem [8]. 

Remarkably, the magnetic excitations at high energies 
are understood less well. But they have come into fo¬ 
cus recently due to remarkable experimental progress in 
resolution [9, 10] and in the direct excitation at high en¬ 
ergies [11]. Conventional spin-wave (magnon) theory in 
powers of 1/5 [12-16] fails to agree quantitatively with 
unbiased quantum Monte Carlo (QMC) simulations [17] 
and high-order series expansions about the Ising limit 
[17, 18]. Even in third order in 1/5, the approach based 
on magnons yields only a weak dip of the dispersion of 
about 3% at (n, 0) compared to its value (7t/2, n/2) while 
series expansion and QMC yield about 10%. Inelastic 
neutron scattering (INS) finds 7(1)% in deuterated cop¬ 
per formate [9, 10]. In analogy to the similar feature in 


the dispersion of excitations in 4 He [19] we call this dip 
‘roton minimum’. 

The agreement of the INS data with the QMC and se¬ 
ries results for the roton minimum on the one hand and 
the strong deviation of spin wave theory on the other 
hand have revived the idea, see, e.g., Ref. 20, that at high 
energies magnons are not the appropriate elementary ex¬ 
citations and that spinons provide a better description 
[10]. The magnons at lower energies would result as con¬ 
fined states of spinon pairs [21] while the spinons are 
deconhned at least at certain points in the Brillouin zone 
(BZ) at higher energies. 

The ambiguity between excitations having fractional 
or integer spin arises in gapless systems because without 
gap it is an ill-posed problem to determine the content 
of quasi-particles. This was for instance seen in one¬ 
dimensional (ID) spin chains [22], But in contrast to 
the 2D case, where magnons are the established Gold- 
stone modes above the long-range ordered ground state, 
the ID Heisenberg chain is known to display fractional 
massless spinons [23]. Consequently, a consistent deter¬ 
mination of the full magnetic excitation spectrum for the 
2D Heisenberg model from low to high energies is still 
missing, but highly desirable. 

This is the central goal of the present Letter. We pro¬ 
vide a quantitative description of the low- and the high- 
energy part of the dispersion in terms of magnons. Tech¬ 
nically, this is achieved by extending continuous unitary 
transformations to continuous similarity transformations 
(CST) which are carried out in momentum space using 
the scaling dimensions of operators as appropriate trun¬ 
cation criterion. 

Our results indicate a strong attractive magnon- 
magnon interaction. This interaction induces a Higgs res¬ 
onance built from two-magnon states, also addressed as 
longitudinal magnon [24-26] . Thereby, significant weight 
in the three-magnon continuum is shifted towards lower 
energies and it can also be seen as magnon-Higgs contin¬ 
uum. These scattering states depress the single magnon 
states and reduce their energy yielding the roton min- 
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imum. In this way, this dip in the dispersion is a fin¬ 
gerprint of strong scattering between magnons and the 
damped Higgs mode. 

2D Heisenberg model — We consider the Heisenberg 
quantum antiferromagnet on the square lattice 

h = jJ2^-s j , ( i ) 

(id) 

where S t denotes the vector operator of a spin 5 = 1/2 
at site i and the sum runs over pairs of nearest neigh¬ 
bors. The ground state displays long-range Neel order. 
Consequently, the classical Neel state with spin up (spin 
down) on sublattice A ( B ) serves as the appropriate ref¬ 
erence state with respect to which bosonic operators cre¬ 
ate deviations. We use the non-hermitian Dyson-Maleev 
representation [27, 28] as is standard in analytic higher 
order calculations [12-16]. Its asset is that the Hamil¬ 
tonian only comprises terms that are at most quartic in 
the bosonic operators. 

In this way, we obtain an initial Hamiltonian first in 
real space and by Fourier transformation in momentum 
space. A Bogoliubov transformation introducing a k ^ and 
/3 k ' ^ as new bosonic operators yields [ 6 , 2! 

H = E 0 + ^ k (a| t a k + /3 k /3 k ) + T + V. (2) 

k 

The momenta k are taken from the magnetic 
BZ containing N points given by the diamond 
with the four edges (± 7 r, 0) and (0, ± 7 r). The 
coefficients E 0 /N = —2 J(S 2 + AS + A 2 /A) and 
cck = 2J(25 + A)^/l — q k with 7 k = (cos k x + cosfc y )/2 
and A := £] k (l — \/l — 7k) ~ 0.158 correspond to 

the ground-state energy per site and the one-magnon 
dispersion in next-to-leading order spin-wave theory 
[ .2, 16]. The Bogoliuov part T := Tq + T° is bilin¬ 
ear; T™ consists of to creation and n annihilation 
operators. The quartic terms are split according to 
V := Vf + Vf + V 3 + Vq + V 4 with the same pattern 
of bosonic operators. By construction, T = 0 in (2) 
once the Bogoliubov transformation is carried out self- 
consistently. We define the part with more (less) creation 
than annihilation operators by H + := Eg + Vf + Vg 

{H~ :=r§ + v! + v2). 

Our goal is to derive an effective Hamiltonian HeS 
which conserves the number of magnons in order to sim¬ 
plify the subsequent analysis. We achieve this by extend¬ 
ing the approach of CUTs [30] to similarity transforma¬ 
tions using the flow equation d{H{€) = [g{tj,T~L {(.)\, where 
£ is a continuous auxiliary variable which labels the stage 
of the CUT: H(£) has the same structure as in (2), but 
all coefficients acquire a dependence on £. Initially, one 
has TL(£ = 0 ) := H and finally H c s '■= H(£ = 00 ). 

We use the generator r/(£) := H + (£) — "H _ (£) [31, 32] 
leading to an effective Hamiltonian H. c s that conserves 


the number of magnons. This is well-established for the 
converging flow of hermitian Hamiltonians and carries 
over to H in (2) because the Dyson-Maleev transfor¬ 
mation does not mix states with different rnagnon num¬ 
ber, but only attributes prefactors to them which spoil 
its unitarity. Summarizing, one maps the complicated 
many-body problem to an effective few-body Hamilto¬ 
nian being block-diagonal in the number of magnons. 
The one-magnon block in particular is diagonal in mo¬ 
mentum space after the CST so that the dispersion J) k 
can be read off. 

Computing [r](£), H(£)] on the right hand side of the 
flow equation leads to a proliferating number of terms. 
One needs a systematically controlled truncation scheme 
to close the flow equation. So far, perturbative trunca¬ 
tions have been used mostly [31, 33-35]. In real space 
approaches the spatial range of processes is a related 
truncation criterion appropriate for gapped systems with 
finite correlation length [32, 35, 36]. In ID, the scaling 
dimension was used in the operator product expansion of 
vertex operators [37, 3< ]. The criterion of the scaling di¬ 
mension is particularly suitable for gapless systems with 
infinite correlation length. Hence, we apply this criterion 
for the 2D massless magnons. 

We consider a generic term T 

T=J + ... + k„)d 2 k!... d 2 k„, (3) 

where C3 k k is a monomial of n bosonic operators of 
creation or annihilation type [39]. If we rescale all mo¬ 
menta k i —> Aki with A < 1 and assume that the coeffi¬ 
cient is a homogeneous function CAki...Ak„ = A c C'k 1 ...k„ 
the term T acquires a prefactor A" -2+c . The scaling di¬ 
mension is the exponent d := n — 2 + c. It quantifies 
the relevance of the term if we zoom to small energies 
in the sense of renormalization. Smaller scaling dimen¬ 
sions stand for more important terms. If the coefficient C 
is not a homogeneous function we determine its dimen¬ 
sional contribution c from the leading term in a series 
expansion in A. 

The bilinear kinetic energy cck(ct: k a: k + /3 k /? k ) has scal¬ 
ing dimension d = 2 — 2 + 1 = 1 because the dispersion 
is linear for small momenta. Thus, this scaling dimen¬ 
sion is the one to compare to; conventionally it is called 
marginal. Terms with smaller dimension would be rel¬ 
evant, terms with larger dimension irrelevant. The Bo¬ 
goliubov terms T = ]T k rk(£)(a k /3L k + h.c.) are marginal 
because we find Tk(£) oc |k| as for the dispersion. The 
quartic terms in V, however, have scaling dimension 2 
because their prefactor is bounded for small momenta, 
i.e., c = 0. Most importantly, the hexatic three-particle 
terms with n = 6 have scaling dimension 4. 

Thus, the dispersion computed in scaling dimension 
1 is the self-consistent mean-field result given below 
Eq. (2). Here we improve this result by performing a 
fully self-similar CST for scaling dimension 2. We make 
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FIG. 1. One-magnon dispersion w(k) in the magnetic BZ. 
Circles depict CST data for L = 16; diamonds, squares, and 
triangles show data from spin-wave theory [15], series expan¬ 
sion and QMC extrapolated to the thermodynamic limit [17], 
Inset: Extrapolation of the dispersion at momenta (tv, 0) and 
(tv/2,tv/2) (circles). The squares depict the series expansion 
and the triangles up the QMC data [17], 


all coefficients of the Hamiltonian (2) ^-dependent, com¬ 
pute the commutators, normal-order the result, and dis¬ 
card the higher hexatic terms to deduce the differential 
flow equations for the coefficients. The flowing Hamil¬ 
tonian, the generator and the differential equations are 
given in the Supplement [29]. The differential equations 
are solved numerically by an adaptive Runge-Kutta al¬ 
gorithm on finite systems with N = L 2 points in the 
magnetic BZ up to L = 16. We stop the flow at values of 
t where the residual off-diagnality, i.e., a measure of the 
norm of the generator [32], has dropped below 10~ 6 J. 
The error due to the finite Runge-Kutta step size is ap¬ 
proximately also 10 -6 J. Furthermore, we checked that 
beyond this I the ground-state energy and the dispersion 
do not change more than 10~ 9 J so that we consider them 
to be converged within good accuracy. 

Low-energy properties — At low energies, the magnons 
are the natural excitations so that we expect that 
our calculation agrees quantitatively with other results, 
e.g., QMC and series expansion [17], which is indeed 
confirmed. We find a ground-state energy per site 
e 0 = E 0 /(2N) = —0.66939(3) agreeing within 5 • 10~ 5 
with e 0 QMC = —0.669437(5) from QMC [40]. Recall 
that the variational Gutzwiller approach yields e[f utz = 
—0.664 corresponding to a gapped state with finite cor¬ 
relation length [10]. 

The data for the spin-wave velocity v agrees as well 
[29]. It has been determined along two different di¬ 
rections in the BZ yielding v\ (along (k,k)) and V 2 
(along (k, 0)). Their linearly extrapolated thermody¬ 
namic values coincide within the error bar yielding v/ J = 


1.675 ± 0.025. This value is in accord with v ( ~ iMG /J = 
1.673 ± 0.007 from QMC [40]. Summarizing, the low- 
energy properties from the CST is quantitatively consis¬ 
tent with the known properties of the Heisenberg model 
on a square lattice. 

High-energy properties — The complete one-magnon 
dispersion for L = 16 is shown in Fig. 1 and compared to 
results from spin-wave theory in order 1/S 3 [15], high- 
order series expansion, and QMC [7]. As discussed 
above, all approaches agree at low energies. 

For the roton minimum at higher energies, this does 
not hold, see Fig. 1. This is seen best in the dispersion 
difference between momenta (7r, 0) and (7r/2,7 t /2). The 
corresponding energies are extrapolated in 1/L in the in¬ 
set of Fig. 1. The roton minimum represents a relative 
dip of (10.5±2.5)% in QMC [17] and (9.5±2)% for series 
expansions [17, 18] while spin-wave theory even in order 
1/S 3 [15] only yields 3.2%. Remarkably, the CST data 
at L = 16 yields a sizable roton minimum of 7%, and the 
extrapolation to L = oo yields 8%. In view of the uncer¬ 
tainties of all approaches, also the high-energy part of the 
dispersion including the roton minimum is quantitatively 
reproduced by the CST in terms of magnons. 

Discussion — As long as there are processes linking 
states with a single elementary excitation, i.e., a quasi¬ 
particle, to states with two or three quasi-particles, the 
quasi-particles decay into continua of two or three quasi¬ 
particles [41]. For collinear Neel order, the Hamiltonian 
allows only for decay into three quasi-particles. Gener¬ 
ally, this hybridization lowers the energy of single quasi¬ 
particle states. A well-studied example is the spin ladder 
with reflection symmetry where the number of triplons 
(the elementary excitations) can only change by an even 
number [42], The triplon dispersion is pressed down by 
the three-triplon continuum. In asymmetric spin ladders 
without reflection symmetry the single triplon states hy¬ 
bridize already with two-triplon continua. Due to the 
larger phase space of two-triplon states this decay chan¬ 
nel has a larger impact [32, 43, 4 ] which is a general 
feature [41]. We conclude that one must understand the 
multi-magnon continuum above the single magnon state 
in order to quantitatively assess the magnon dispersion 
in general and the roton minimum in particular. 

In gapless systems such as the 2D Heisenberg model 
(1) the continua start just above the single-magnon ener¬ 
gies because multi-magnon scattering states can be built 
from the single magnon and one or two magnons arbi¬ 
trarily close to u >o = 0. Thus, the crucial impact of the 
continua nearby in energy is mainly driven by the low 
energy physics defined by terms of low scaling dimen¬ 
sion. Therefore, the scaling dimension is an appropriate 
criterion even for a calculation concerned with the high- 
energy roton minimum. 

The energy lowering due to the hybridization with 
three-magnon states is less effective as argued in the 
analogous system of spin ladders. But an attractive in- 
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teraction among the magnons shifts spectral weight to 
lower energies and enhances the impact of the magnon 
decay. The marked roton minimum indicates that this 
is a very important aspect. Indeed, the term Vf con¬ 
tains the attractive interaction fH 4 ) between the a and 
P magnons living on the A (S/ ot = —1) and the B sub¬ 
lattice (5* ot = +1), respectively [29]. 



FIG. 2. Two-magnon part of the longitudinal dynamic struc¬ 
ture factor S zz ( k, w). Dotted black line: without interaction, 
dashed blue line: bare interaction, solid red line: renormal¬ 
ized interaction. The curves are computed for L = 192 using 
the matrix elements of the L = 16 system and broadened by 
0.013J. The dispersion is interpolated smoothly, the interac¬ 
tion in a piecewise constant way. 

The crucial effects of the magnon-magnon attraction 
are depicted in Fig. 2 where we computed the longi¬ 
tudinal dynamic structure factor S zz (k, uj) in the two- 
magnon S/ ot = 0 channel without and with (renormal¬ 
ized) interaction. The CST-renormalization of the ob¬ 
servable is not included to focus on the pronounced inter¬ 
action effects. The attraction leads to a shift of spectral 
weight to lower energies and, strikingly, to the formation 
of a resonance which we identify as the Higgs resonance. 
The position of the resonance at ( 7 r, 0) is found to be 
lower than the one at (n/2, 7t/2). This is in line with the 
argument that the single-magnon dispersion is pushed to 
lower energies more strongly at (7T,0) than at (7t/2,7t/2) 
leading to the roton minimum because due to the Higgs 
resonance the decay of a single magnon not only com¬ 
prises the decay into three rather independent magnons, 
but also the decay into two constituents: magnon plus 
Higgs resonance. In accord with previous studies [32, 41- 

] this is an efficient mechanism to lower the magnon 
dispersion and to induce the roton minimum. 

The renormalization of the interaction in the CST flow 


enhances the effect of the attraction further. Though this 
effect is small in absolute terms it makes itself felt in the 
size of the roton dip. If we switch off the flow of the parti¬ 
cle conserving interactions, i.e., without renormalization 
of the magnon-magnon interaction, the dip of the roton 
minimum is reduced to 5%. This also explains why the 
third order calculation in 1/S yielded only a small dip 
although some interaction effect was included (diagram 
(k) in Ref. [15]). Still, the roton minimum is insufficiently 
captured. Note that the renormalization of the interac¬ 
tion can be as large as 50% for certain momenta [29]. 

We emphasize that our line shapes in the longitudi¬ 
nal channel agree remarkably well with the experimental 
data, see Figs. 2c and 2g in Ref. [10]. This holds for 
the position of the Higgs resonances and for the rela¬ 
tive heights. Investigating the experimental transverse 
dynamic structure factor (Figs. 2b and 2f in Ref. [10]) 
one finds pronounced peaks at (7r,0) and at (7 t/2 ,7 t/2 ). 
We attribute both of them to magnons and the smaller 
weight at (n, 0), combined with a continuum tail, to 
strong magnon-Higgs scattering. This also shifts the dis¬ 
persion to lower energies at (71% 0). Hence, the roton min¬ 
imum in the magnon dispersion is a fingerprint of the 
Higgs mode and of strong magnon-Higgs scattering. 

Summary — For the S = 1/2 Heisenberg antiferro- 
magnet on the square lattice, this letter provides a fully 
consistent and quantitative picture in terms of magnons 
which describes the low- and the high-energy physics on 
equal footing. It is achieved by extending the approach 
of continous unitary transformations to continuous simi¬ 
larity transformations. In momentum space, the flow of 
the renormalized couplings is closed by the truncation at 
the level of terms with scaling dimension 2. 

The striking agreement of our findings with the exper¬ 
imental results [10] underlines the validity of the magnon 
quasiparticle picture for long-range ordered quantum 
magnets at all energies. In this way, the importance of 
a significant attractive interaction between magnons of 
different S z ot is established. This interaction induces a 
resonance corresponding to the longitudinal magnon or 
Higgs resonance of the symmetry broken phase. The de¬ 
cay of a magnon into three magnons is strongly enhanced 
because it is effectively a decay into a magnon and a 
Higgs resonance. The roton minimum is the fingerprint 
of this mechanism. 

Further studies should address other response func¬ 
tions providing predictions to spectroscopic experiments. 
We highlight that the developed approach is applicable 
to a humongous variety of long-range ordered phases. It 
will help to identify their dynamical properties and even¬ 
tually their instabilities and ensuing transitions towards 
other exotic quantum phases. 
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SUPPLEMENTARY MATERIAL 


This supplementary materials gives all specific informations concerning the extrapolation of the spin wave velocity 
and the renormalization of the interaction. In addition, details of the CST performed in the main body of the 
manuscript are given. 



6 



t 


FIG. 3. Extrapolation of the results in 1 jL and comparison to other findings. All lines are linear fits. Upper panel. Spin-wave 
velocity v deduced from two different directions V\ (triangles left) and V2 (triangles right) shown in the magnetic Brillioun 
zone (inset). Lower panel Renormalization of the attractive interaction of two magnons for smallest, non-zero k values along 
direction ui (see inset upper panel) vs. the flow parameter i\ results for direction V 2 are almost the same. 


Spin wave velocity and magnon-magnon interaction 

Our data for the spin-wave velocity v agrees well with previous results as shown in the upper panel of Fig. 3. It 
has been determined along two different directions in the BZ yielding v\ and ^ 2 , see upper panel of Fig. 3. Their 
linearly extrapolated thermodynamic values coincide within the error bar yielding v/J= 1.675 ± 0.025. This value is 
in accord with n^ MC /J = 1.673 ± 0.007 from QMC [40]. 

In the lower panel, we illustrate that the CST procedure renormalizes the interaction sizably. For the smallest 
non-vanishing momenta the interaction is attractive and renormalized by up to 50%. Averaged over all interaction 
coefficients the renormalization is of the order of 10%. Note in addition that not all terms of the interaction act 
attractively. 


Flow equation of the continuous similarity transformation 

Below, we give all specific informations concerning the CST performed in the main body of the manuscript: the 
initial Hamiltonian 'H and the initial generator 77 , the flowing Hamiltonian %{€) and the flowing generator r)(£) as well 
as the corresponding flow equations. 


Initial Hamiltonian and initial generator 

The initial Hamiltonian before the CST is given by 

U = E 0 + w k(aka k + /? k /3 k ) + V ■ 

k 


( 4 ) 
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The coefficients of the quadratic part are E 0 /N = — 2 J(S 2 + AS + A 2 /A), Wk = 2J(2S' + A)^/l — with 
7k = (cos k x + cosky) /2, A:= ^D k (l ~ \A — 7k)- The two-magnon interaction V is given explicitly by 


V = -^^<5 g (1 + 2- 3-4) hi 2 l 3 l 4 


N 


1234 


^1234*^1*^2^3^4 "b 2V r 1 ^234Q'i/?_2 a 3Q ; 4 + ~b 4 /^_2 


+ 2V 1 ( 2 5 ; ] 4 /3l 4 a3/3-2/3_i + 2V^_^_ 3 a\P-i + + V^_^_ 2 a 3 a A + .(5) 


The subscripts i = 1, 2, 3,4 stand for the momenta kj and —i stands for — kj. The conservation of momentum in the 
lattice is ensured by the Kronecker symbol Sq (k) which implies k = 0 modulo reciprocal lattice vectors from the 
reciprocal lattice T^ of the H-sites, i.e., g € T^ means g = ( mr,rmr ) with the integers n,m if the lattice constant of 
the original square lattice is set to unity. The vertex functions Vy 234 are given explicitly in Ref. 16. 


Flowing Hamiltonian and flowing generator 


We include all terms in the flowing Hamiltonian H(i) with a scaling dimension d < 2. Explicitly, one finds 
7~L{1) = E 0 + w (1) {a\a 1 + + T (1) 


+ 


El, 2 , 3,4 {Cl (1234) a{a\a 3 a A + C 2 (1234) a[a 2 ^ 4 + C 3 (1234) 

C 4 (1234) a\a\a z p\ + C 5 (1234) 4 + C 6 (1234) a[a 2 a 3 /. 3 4 


+ 


+ C 7 (1234) + C 8 (1234) + C 9 (1234) a^P^} 

The corresponding flowing generator is then given by 

V(t) = Ei r (!) - “lEi) 

+ El, 2,3,4 {C4 (1234) a[aia 3 pi + C 5 (1234) a^sAi ~ <? 6 (1234) a{a 2 a 3 p 4 
- C 7 (1234) ai 4p 3 /3 4 + C 8 (1234) aW 2 PM~C 9 (1234) ai a 2 /3 3 P 4 J . 


( 6 ) 


( 7 ) 


The coefficients E 0 , oj (1), T (1), and Ci (1234) depend on the flow parameter t and satisfy the initial conditions 


E 0 \ e=Q = -2J(S 2 + AS + A 2 /4) 

W (1)|, =0 = 2J(25 + H)^1- 7i 2 

r ( 1 )L=o = 0 

Ci (1234) | fcQ = -hhhU^V^6 G (1 + 2 - 3 - 4) 

C 2 (1234) | fcQ = -^hhhh^V^i^G (1 - 2 + 3 - 4) 

C 3 (1234) | f=Q = -hhhh^V^_ m 8 G (1 + 2 - 3 - 4) 
C 4 (1234) | /=Q = -‘IhhhU^V^c (1 + 2 - 3 + 4) 

C 5 (1234) | fcQ = -^khh^vf^G (1 + 2 + 3 - 4) 
C 6 (1234) | f=Q = -2 hhhU^V^c (1 - 2 - 3 - 4) 

C 7 (1234) | f=Q = -2hhhl 4 ^vf 31 . 2 6 G (-1 + 2 - 3 - 4) 
C 8 (1234) |^_ Q = — hhhhj^ ^ 2 - 3 - 4 ^° (1 + 2 + 3 + 4) 

C 9 (1234) |, =0 = -hhl 3 l 4 ^V_f A12 S G (-1 - 2 - 3 - 4) 


(8a) 

(8b) 

(8c) 

(8d) 

(8e) 

(8f) 

(8g) 

(8h) 

(8i) 

(8j) 

(8k) 

( 81 ) 







Note that the hexatic three-magnon part of the Hamiltonian consisting of six magnon annihilation or creation 
operators has a scaling dimension d = 4 and are therefore not relevant at studied level of truncation d = 2. 


Flow equations 


Inserting 'H(l) and r/(£) into the flow equation dfT-L(£) = 'H(i)] and keeping self-similarly all operators already 

present in ?£(£), one obtains the following flow equations 


dlEo = El, 23, 4 (-8)Cs (1234) C 9 (1234) + £(-2)r (1) T (1) (9a) 

1 

cW(l) = (-2)r(l)r(l)+ E 3 ,4 {-4<7 4 (3415)<7 6 (1345)<5 g (3 + 4-1 + 5) 

-16(78 (1345) <7 9 (1345) S G (3 + 4 + 1 + 5) 

-4<7 8 (131-3) T (3) - 4<7 6 (113-3) T (3)| (9b) 

8iT (1) = -2r (1) u (1) - 8 E 3 ,4 >5 {<^4 (3415) Cg (34-15) <5 G (-1 + 3 + 4 + 5) 

—8(7 5 (345-1) C 9 (1345) S G (1 + 3 + 4 + 5) 

E,3 {(72 (13-1-3) T (3) 

(7 8 (13-1-3) T (3)| (9c) 

diCi (1234) = <5 G (1 + 2 — 3 — 4) {(—1)(7 4 (1243) T (3) 

(—1)<7 4 (1234) r (4) 

(-1)0, (2345) T(l) 

(—1 )<7 6 (1345) r (2) 

Es, 6 (-4)^4 (2546) C 6 (1356) S G (2 + 5 — 4 + 6) 

(— 4)<7 4 (1536) C 6 (2456) 6 G (1 + 5 - 3 + 6) 

(—4)<7 8 (1256) C g (3456) S G (1 + 2 + 5 + 6) j 

(10a) 


diC 2 (1234) = <5(3 (1 — 2 + 3 — 4) {(—1)(7 4 (1-423) T (-4) 

(—1)<7 5 (13-24) r (-2) 

( —1)<7 6 (12-34) r (-3) 

(-I)CV (234-1) r(-l) 

E 5 ,6 (-4)^4 (5623) Ce (1564) <5 G (5 + 6 - 2 + 3) 

(—4)<7 5 (1564) C 7 (2356) <5 G (1 + 5 + 6 - 4) 

(—8)<7 4 (1526) C 6 (5346) <5 G (1 + 5 - 2 + 6) 

(—8)<7 5 (5364) C 7 (1256) <5 G (5 + 3 + 6 - 4) 

(—32)(7 8 (1536) Cg (2546) S G (1 + 5 + 3 + 6) j 


(10b) 
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9,(73(1234) = J G (l + 2-3-4){(-l)C 5 (-3124)r(-3) 

(-1)C 5 (-4123) T (-4) 

(—1)CV (-2134) r (-2) 

(-1)^7 (-1234) T(-l) 

Es (-4 )C 5 (5264) C 7 (5136) S G (5 + 1 + 6 - 4) 

(—4)<7 5 (5163) C 7 (5246) S G (5 + 1 + 6 - 3) 

(—4)<7 8 (5612) C 9 (5634) S G (5 + 6 + 1 + 2) j (10c) 


9,(74 (1234) = (5 G (1 + 2 — 3 + 4) j(w (3) — u> (1) — u (2) — ui (4)) C 4 (1234) 

(—2)<7i (123-4) T (-4) 

(C 2 (234-1) r ( 1 ) + C 2 (134-2) T ( 2 )) 

(—4)(7 8 (124-3) r (3) 

E 5 , 6 (—2)<7i (1256) <7 4 (5634) S G (1 + 2 - 5 - 6) 

(—1)<7 2 (2546) Ci (1536) S G (2 - 5 + 4 - 6) 

(—1)<7 2 (1546) (7 4 (2536) <5 G (1 - 5 + 4 - 6) 

(—8)<7 6 (2356) (7 8 (1546) <5 G (2 - 3 - 5 - 6) 

(—8)<7 6 (1356) C 8 (2546) 5 G (1 - 3 - 5 - 6) 

(—4 )C 7 (3456) <7 8 (1256) <5 G (3 - 4 + 5 + 6)} (lOd) 


9,(7 5 (1234) = $ G (1 + 2 + 3 — 4) j(—w (1) — w (2) — uj (3) + u> (4)) C 8 (1234) 

(-2)<7 3 (234-1) T (1) 

(<7 2 (1-234) T (-2) + C 2 (1-324) T (-3)) 

(-4)<7 8 (1-423) T (-4) 

Es,6 (— 2)<7i (5623) C 6 (1564) 6 G (5 + 6 - 2 - 3) 

(—1)<7 2 (1536) C 5 (5264) 5 G (1 - 5 + 3 - 6) 

(—1)<7 2 (1526) C 5 (5364) S G (1 - 5 + 2 - 6) 

(—8)CV (5346) C 8 (1526) 6 G (5 - 3 + 4 + 6) 

(—8)CV (5246) C 8 (1536) 6 G (5 - 2 + 4 + 6) 

(-4)<7 6 (1564) C 8 (5623) 5 G (1 - 5 - 6 - 4) j (lOe) 
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diC e ( 1234) = (5 G (1 — 2 — 3 — 4) |(w (1) — ui (2) — u> (3) — u> (4)) C” 6 (1234) 

(—2)Oi (1-423) T (-4) 

(0 2 (13-24) T (2) + C 2 (12-24) T (3)) 

(-4)C 9 (234-1) T(l) 

E 5 ,e (—2 )Ci (5623) O e (1564) 4 G (5 + 6 - 2 - 3) 

(-1)0 2 (5364) O e (1256) <5 G (5 - 3 + 6 - 4) 

(-1)0 2 (5264) O e (1356) S G (5 - 2 + 6 - 4) 

(—8)C 4 (1536) Og (2546) 6 G (1 + 5 - 3 + 6) 

(-8)04 (1526) O g (3546) 4 G (1 + 5 - 2 + 6) 

(— 4)05 (1564) Og (2356) 4 G (l + 5 + 6-4)| (lOf) 

diC 7 { 1234) = <5 G (1 — 2 + 3 + 4) |(—w (1) + w (2) — w (3) — w (4)) O 7 (1234) 

(-2)0 3 (2-134) T(l) 

(C 2 (-3124) T (-3) + 0 2 (-4123) T (-4)) 

(—4)Og (1-234) r (-2) 

Es,6 (—2)0 3 (5634) 0 7 (1256) S G (5 + 6 - 3 - 4) + 

(-1)0 2 (5164) O 7 (5236) S G (5 - 1 + 6 - 4) 

(-1)0 2 (5163) O 7 (5246) 6 G (5 - 1 + 6 - 3) 

(—8)0 5 (5264) Og (1536) S G (5 + 2 + 6 - 4) 

(—8)0 5 (5263) Og (1546) S G (5 + 2 + 6 - 3) 

(—4)0 4 (5612) Og (5634) S G (5 + 6 - 1 + 2 ) j (10g) 

diC s ( 1234) = <5 G (1 + 2 + 3 + 4) |(—w ( 1 ) — oj ( 2 ) — oj (3) — ui (4)) Og (1234) 

Es ,6 (— 2 )Oi (1256) O s (5634) S G (1 + 2 - 5 - 6 ) 

(-1)0 2 (2546) O s (1536) S G (2 - 5 + 4 - 6 ) 

(-1)0 2 (1546) O s (2536) (5 G (1 - 5 + 4 - 6 ) 

(-1)0 2 (2536) O s (1546) <5 G (2 - 5 + 3 - 6 ) 

(-1)0 2 (1536) O s (2546) 5 G (1 - 5 + 3 - 6 ) 

(—2)0 3 (3456) O s (1256) ,5 G (3 + 4 - 5 - 6 )| (lOh) 

a ; O g (1234) = <5 G (1 + 2 + 3 + 4) |(—w (1) — oj ( 2 ) — uj (3) — oj (4)) Og (1234) 

E 5 ,6 (— 2 )(— 2 )Oi (5612) Og (5634) 5 G (5 + 6 - 1 - 2 ) 

(-1)0 2 (5264) Og (1536) S G (5 - 2 + 6 - 4) 

(-1)0 2 (5164) Og (2536) S G (5 - 1 + 6 - 4) 

(-1)0 2 (5263) Og (1546) ,5 G (5 - 2 + 6 - 3) 

(-1)0 2 (5163) Og (2546) 6 G (5 - 1 + 6 - 3) 

(—2)0 3 (5634) Og (1256) S G (3 + 4 - 5 - 6 ) | (lOi) 

Note that the Kronecker symbols S G (k) in the sums are redundant since the coefficients intrinsically conserve the 
total momentum by definition (see Eq. ( 8 )). They are included to underline momentum conservation. 


